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We present an overview of our understanding of electrostatic interactions between charged macro- 
(3J[), molecular surfaces mediated by mobile counter- and coions. The dichotomy between the weak and 

the strong coupling regimes is described in detail and the way they engender repulsive and attractive 
' interactions between nominally equally charged macroions. We also introduce the concept of dressed 

counterions in the case of many-component Coulomb fluids that are partially weakly and partially 
strongly coupled to local electrostatic fields leading to non-monotonic interactions between equally 
charged macroions. The effect of quenched surface charge disorder on the counterion-mediated 
electrostatic interactions is analyzed within the same conceptual framework and shown to lead to 
fZt unexpected and extraordinary electrostatic interactions between randomly charged surfaces with 

^ ' equal mean surface charge densities or even between effectively neutral macroion surfaces. As a 

result, these recent developments challenge some cherished notions of pop culture. 
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I. INTRODUCTION 



The nature of electrostatic interactions between charges started to be studied intensively in the second half of the 
18th century Benjamin Franklin first inferred from an observation that surprisingly there is no force on a charge 
inside a charged sphere, an observation later repeated by Joseph Priestley (1767). John Robison (1769) determined 
that the electrostatic force falls off with (almost) the second power of separation between charges. Based on the work 
of Charles A. de Coulomb (1777) who is a co- inventor of the torsion balance for measuring the force of magnetic 
and electrical attraction, Henry Cavendish measured directly the interactions between charges (1779) but did not 
{y-^ | publish his results. They were eventually published by William Thomson (Lord Kelvin) one hundred years after the 
£f) . original discovery (1879). From hard science electrostatic interaction penetrated the pop culture in general and today 

' everybody knows that opposites attract and likes repel! 
qq , The exact form of the electrostatic interaction is since known to be given by the Coulomb's law which states that 
the interaction potential between two charges e\ and ei in vacuo located at r and r' respectively, can be written in 
^> \ the standard form (in SI units) as 



were eo is the permittivity of vacuum. Electrostatic interaction is the fundamental interaction in molecular world 
giving rise to short range intramolecular bonds as well as longer ranged interactions between molecules and their 
aggregates Q. While atomic bonds are obtained by combining the Coulomb potential with the quantum of action, 
longer ranged colloidal and nanoscale electrostatic interactions are obtained by combining the Coulomb potential with 
the thermal energy k^T . 

Almost exactly a hundred years ago Gouy Q and Chapman Q were the first to combine thermal energy and 
Coulomb interactions into a statistical theory of Coulomb fluids basing their approach on what became later known 
as the Poisson-Boltzmann (PB) equation. Coulomb fluid in general is an assembly of (variously) charged particles 
in thermal equilibrium. This mean-field approach to Coulomb fluids was gradually elaborated in uncountably many 
ways starting from the seminal work of Debye and Hiickel [j| and finally codified as a cornerstone of the fundamental 
DLVO theory of colloidal interactions by Derjaguin and Landau as well as Verwey and Overbeek 0], where it 
figures as the repulsive electrostatic part of the total disjoining pressure acting between macromolecular bodies [1, [§] . 
The attractive component in this case is provided solely by the van der Waals interactions that have their origin in 
the quantum and thermal fluctuations of electromagnetic fields [Toj ] . 

One could claim that modern formulation of statistical mechanics of Coulomb fluids starts with the work of Edwards 
and Lenard [Til ] where the grand-canonical partition function of a Coulomb fluid was written in the form of a functional 



2 



integral over fluctuating electrostatic fields. While Podgornik and Zeks [l2| realized that the collective description 
based on the mean-field PB theory arises from the saddle-point approximation to this functional integral in the case 
of a weak coupling regime corresponding, for instance, to small external charges, Netz showed (l3l . Il4 | that another 
approximation valid in the regime of large external charges leads to a completely different fixed point, formulated in 
terms of a single-particle description of the same system. This so-called strong coupling regime has been the focus 
of various theoretical studies over the past decade [l5l - [T9l |. The two approximations were shown to correspond to 
extremal values of a single electrostatic coupling parameter and to consistently bracket all available simulation results 
on statistical properties of non- homogeneous Coulomb fluids flM [l6j. i.e., Coulomb fluids confined between charged 
boundaries. 

In what follows we shall describe the main consequences of this weak-strong coupling dichotomy especially as 
they transpire in the case of interactions between charged (planar) macromolecular surfaces across a Coulomb fluid 
comprising mobile counter- and coions. We shall see that the counterion-mediated electrostatic interactions are strictly 
repulsive in symmetric external charge configurations for the case of weak coupling but can turn attractive for strongly 
coupled surfaces. In the case of many-component systems, composed typically of a weakly coupled univalent salt and 
strongly coupled polyvalent counterions, these interactions show a subtle interplay between repulsions and attractions 
for nominally equally charged surfaces. On the other hand, the presence of quenched charge disorder on bounding 
surfaces can lead to pronounced electrostatic attraction. These exotic features only arise in special cases of strongly 
coupled systems and can not be rationalized within the standard mean-field PB approach. 

II. SCENERY 

Electric charges and electrostatic interactions are ubiquitous in soft- matter and biological systems [2(| Hl[ . Soft 
materials are typically composed of macromolecules such as polymers, colloids and proteins which often acquire surface 
charges when dissolved in a polar solvent like water. This is usually due to dissociation of surface chemical groups, 
which leaves permanent charges on macromolecular surfaces and releases oppositely charged microscopic counterions 
into the solution. Soft materials are easily deformed or re-arranged by interaction potentials comparable in magnitude 
to thermal energy. It thus becomes clear that electrostatic interactions, that are typically long-ranged and strong, 
constitute a prominent factor in determining the behavior and properties of soft materials. This makes charged 
materials central to many technological applications and on the other hand, a challenging subject for fundamental 
research in inter-disciplinary sciences 0. In what follows, we briefly review a few examples to demonstrate the 
diversity of phenomena associated with charged soft-matter systems. 

A. Colloids, polymers and membranes: The mesoscopic scale 

Colloids are abundant in nature and industry: smoke, fog, milk, paint and ink are only a few examples of colloidal 
systems. They comprise tiny solid or liquid particles that are suspended in another medium such as air or another liq- 
uid. An important factor, which makes colloidal solutions in many ways different from molecular or simple electrolyte 
solutions (such as sugar or salt solution), is the large asymmetry in size and mass between the colloidal particles and 
solvent molecules (or microscopic ions): colloids are mesoscopic or even nanoscopic particles with sizes in the range 
of a few nanometers to microns that are indeed made of many atoms, but not yet sufficiently many to make them 
behave like macroscopic bodies. 

In colloidal dispersions, the total area that is in contact with solvent is tremendously large: for nanometer-sized 
colloids, nearly half of the atoms are at the surface whereas for macroscopic bodies this ratio tends to zero. This is even 
more true for extended quasi-two-dimensional macromolecular aggregates such as lipid membranes and surfactant- 
covered interfaces that can carry a substantial amount of charged molecular moieties [2l| . Therefore, contrary to the 
typical situation in the macroscopic world, the physics of mcso- and nanoscopic particles are dominated by "surface" 
properties and interactions @-[9|. 

Another relevant mesoscopic or macromolecular system are polymers (with everyday- life examples like chewing gum, 
dough or egg white), in which many repeating subunits (monomers) are chemically connected to form a flexible chain. 
Flexible polymers are distinguished by their many degrees of freedom associated with conformational rearrangements 
of monomers that are easily excited by thermal energy at room temperature leading to a diverse phase behavior 
spanning extended and strongly entangled polymers such as DNA in solution, all the way to collapsed polymer 
chains organizing into compact globular states as in the case of proteins. Depending on their chemical structure, 
polymer chains can have a large mechanical stiffness as well, behaving like rigid rods at small length scales, or can 
be substantially charged giving rise to poly electrolytes, in both cases playing an important role in biological processes 
occurring in the cell |22| . 
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FIG. 1: Schematic view of a charged solution consisting of charged colloids (dark spheres), charged polymer chains and 
microscopic neutralizing counterions (small black spheres). Colloids may be covered by charged polymer brushes, which 
generate an additional repulsive interaction between them. 



B. Charges: from industry to biology 

In the mesoscopic world, only electromagnetic interactions that emerge in a variety of forms, are important; they 
are capable of overcoming thermal fluctuations-which are characterized by an energy scale of about 0.025eV at room 
temperature 0-and thus enable formation of stable condensed phases for soft complex materials. 

In general, colloids dissolved in an aqueous solvent attract each other due to ubiquitous van der Waals dispersion 
forces that result from induced charges on their surface @- As a result, colloidal particles tend to form large aggregates 
that typically sediment and destroy the dispersion. In many applications (for example in food emulsions such as milk), 
however, stability of a colloidal dispersion is a desirable property. One way to stabilize dispersions against aggregation 
is to generate long-range repulsive interactions between colloidal particles by charging their surfaces, which leads to 
the DLVO mechanism for the stability of colloidal dispersions g, 0] ■ Another method of stabilization is to end-graft 
polymer chains (or polymer brushes) to the particle surfaces [23j . For this task, charged polymers are ideal since they 
swell substantially in aqueous solutions and inhibit close contact between colloids (see Fi g, [p . This latter mechanism 
has the advantage that it is highly stable against the addition of electrolyte or salt ions |24| . 

Charged polymers, or poly 'electrolytes, play a significant role in the production of cheap, non-toxic and environmen- 
tally friendly materials (HI, H^]. In contrast to water- insoluble hydrocarbon chains, polyelectrolytes typically show 
high affinity for water and heavy metal ions, which makes them useful in applications such as super-absorbing diapers, 
waste water purifiers and washing detergents and their additives. 

In biology electrostatic effects between charged polymers such as DNA and RNA are ubiquitous |2(j. DNA, for 
instance, is a long biomolecule with a total length of about two meters in human cells, bearing one elementary charge 
per 1.7A, which for human DNA adds up to 10 10 elementary charges overall! Yet the DNA is densely packed inside the 
cell nucleus with a diameter of less than a few microns. In cucaryotic cells, this storage process involves a hierarchical 
structure on the lowest level of which short segments of DNA are tightly wrapped around positively charged histone 
protein complexes of a few nanometers in diameter [22| . This protein-DNA complexation is believed to be governed 
by electrostatic interactions [HI, H3, HH . Electrostatic effects also play a key role in complexes of DNA with cationic 
lipids [29M3~H , which are promising synthetically based non- viral carriers of DNA for gene therapy [Hf • 

Another example (which is closely related to the results presented later in this chapter) is the DNA condensation 
[20I [33I HH , in which electrostatic effects enter in a counter- intuitive fashion: here like-charged segments of DNA 
strongly attract each other! In the in vitro experiments [35j . the condensation of DNA is realized using bacteriophages, 
which consist of a rigid shell (the capsid) that accommodates a single molecule of viral DNA. These viruses can inject 
their DNA into a cell or a lipid vesicle. As a result, large lengths of DNA (up to a hundred microns) can be fitted and 
condensed into a tightly packed, circumferentially wound torus with a diameter of about a hundred nanometers. This 
packaging process, which works against the Coulomb self-repulsion and the conformational entropy of the DNA chain, 
is facilitated and depends upon the presence of high- valency counterions in the medium. Similarly, other highly charged 
polymers, such as negatively charged F-actin and microtubules can aggregate into closely packed rod-like bundles when 
small amounts of polyvalent cations are added to the solution |36l - l38j . It turns out that, in general, when particles 
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are strongly charged, the role of electrostatic interactions dramatically changes (l5j: here electrostatic interactions 
themselves can trigger the destabilization of charged solutions by mediating attractive like-charge interactions! 



C. Theoretical challenge and coarse-grained models 



From a theoretical point of view, charged systems pose a many-body problem: macroions, such as charged colloids 
and polymers, and other charged macromolecular surfaces, such as lipid membranes and surfactant layers, are always 
surrounded by countcrions, and also in general by coions. These particles form loosely bound ionic clouds around 
macroions and tend to screen their charges. In particular, countcrions that are attracted towards macromolecular 
surfaces, predominantly determine the static and dynamic properties of macromolecular solutions. Understanding 
the interactions between macromolecules across an ionic medium thus requires an understanding of the countcrionic 
clouds first. 

In the most common theoretical approaches known also as primitive models, the molecular nature of the solvent is 
neglected and is represented by a continuum dielectric medium. In reality, the solvent structure is locally perturbed 
around particles, which can give rise to additional short-ranged solvent-induced interactions 0, H^, H(|. On the other 
hand, the microscopic features of the macroions are taken into account using coarse-grained models that incorporate 
only a few effective parameters such as an effective surface charge density. In most cases, the specific effects associated 
with ions [1^] as well as the image charge effects due to dielectric inhomogeneities are also neglected. These models 
therefore represent a crude simplification of reality, yet given those simplifications, they can still lead to a systematic 
and clear understanding of electrostatic effects. 

Here we shall first begin by adopting such a simple model for the interaction between charged macromolecular 



surfaces in the presence of counterions only, but then examine the effects due to the additional salt [4 11 ] and the 
heterogeneous or disordered distribution of surface charges [43 - |45j in more detail. For simplicity, we shall also focus 
only on the case of planar charged surfaces, appropriate for the case of char ged membranes, solid substrates, or large 
colloids. Other factors such as non-planar geometry of charged surfaces [46N48j . image charges [43l. |48M50| . dissimilar 
surfaces [5lj ] or multipolar structure of counterions (52j have been analyzed within the same context. 



III. LENGTH SCALES IN A CLASSICAL CHARGED SYSTEM 



Consider a system of fixed charged objects with uniform surface charge density — (7 s en (with eo being the elementary 
charge) that are surrounded by their neutralizing countcrions of charge valency +q in a solvent of dielectric constant 
e at temperature T. 1 

One of the basic characteristic length scales in a charged system is the Bjerrum length [53[ 

£ B = e 2 J(4iree k B T), (2) 

which is set by the ratio between the thermal energy, k B T , and the Coulomb interaction energy between two elementary 
charges at separation r, i.e., V/(k B T) = t^/r. The Bjerrum length thus measures the distance at which two elementary 
charges interact with an energy equal to k B T. In water and at room temperature (e = 80), one has £ B — 7.lA. For 
counterions of charge valency +q, the Bjerrum length may be redefined as q 2 l B . 

Other length scales may be identified by considering the specific form of the charge distribution and geometry of 
macroions. For uniformly charged planar surfaces (Fig. [5]), one can define another key length scale by comparing the 
thermal energy with the energy scale of the counterion-wall attraction, i.e., uj ' (k B T) = z/fi, where z is the distance 
from the wall and 



/! = l/(2irq£ B a s ) (3) 

is known as the Gouy- Chapman ( GC) length [1,0]. The GC length measures the distance at which the thermal energy 
equals the counterion-wall electrostatic interaction energy. It also gives a measure of the thickness of the counterion 
layer at a charged wall as we shall see later. 



We conventionally assume that macroions are negatively charged and counterions are positively charged, thus cr B and q are both positive 
by definition. 
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FIG. 2: Schematic view of the structure of a layer of counterions at an oppositely charged surface for small (a) and large (b) 
value of the coupling parameter. 



For planar systems where no other length scales are present, it follows that only the dimensionless ratio between 
the above two length scales matters, i.e., 



5 = q 2 £ B /p = 2Trq 3 £la a . 
This parameter is known as the electrostatic coupling parameter [l3l . [l4| . 



(4) 



IV. FROM MEAN-FIELD TO STRONG COUPLING REGIME 



For small coupling parameters 3 <C 1, equation (|4|) shows that the GC length is relatively large, which indicates 
that counterions form a loosely bound cloud at an oppositely charged wall (Fig. For large coupling parameter 

E 3> 1, in contrast, the GC length is relatively small and counterions are strongly attracted toward the wall (Fig. 
[2}d). Further insight may be obtained by considering the typical distance between counterions at a charged surface. 
For counterions residing near the surface, the local electroneutrality condition implies a typical lateral separation of 



since each counterion neutralizes the charge of an area that scales as a\ ~ q/a s 
independent length scale and may be written in the units of the GC length as 



(5) 

Counterion spacing aj_ is not an 



(6) 



A. Weak coupling or mean-field regime 



In the regime where S< 1, equation (|6|) shows that the lateral separation of counterions near the surface is small 
compared with the typical layer thickness, p, which further indicates that the counterions tend to form a diffuse fluid- 
like layer at the surface (Fig. [2k). 2 This regime is dominated by collective mean- field-like effects, i.e., counterions 
become uncorrelatcd from each other in a statistical sense as each counterion in the layer interacts with a diffuse 
cloud of many other counterions. Therefore, S< 1 identifies the weak coupling (WC) or mean-field regime, which is 
relevant to systems with weakly charged surfaces, low valency counterions and/or high temperature p~5| ] . 

Formally, one can employ a mean-field approximation in order to describe the system in the WC regime by neglecting 



all inter-particle correlations on the leading order. The mean-field approximation is exact in the limit S — > 54j and 



leads to the so-called Poisson-Boltzmann (PB) equation for the mean-field electrostatic potential ip(r) 0, H[, i.e. 



££o 



p (r) + qe n Q(r) cxp(- Pqe i/)) , 



(7) 



where @ = l/(k B T). 
nal charges (macroions) 



The first term on the r.h.s. represents the charge distribution due to fixed exter- 
Po{t), and the second term is the mean- field number density of counterions npB(r) = 



A more accurate estimate of the typical distance, a, between counterions in an extended three-dimensional layer gives a/fi ~ r. 1 / 3 |l3t . 
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TABLE I: Typical values of physical parameters for some realistic systems. 



uqVI(y) exp(—f3qeot/j), where no is a normalization prcfactor. The "blip" function f2(r) is equal to one in the region 
accessible to countcrions and zero elsewhere. 

For point-like counterions at a single uniformly charged wall, the PB equation yields the well-known algebraically 
decaying density profile 0, Q 

^pb(z) _ 1 n 

where z is the distance from the wall. Note that the density profile is normalizable to the total number of counterions 
in order to ensure electroneutrality, and that the GC length, a, equals the height of a layer containing half of the 
counterions, thus giving a measure of the typical layer thickness. The contact density of counterions npe(z = 0) = 
27r^BCs turns out to be an exact result within the present model and remains valid beyond the mean- field theory [55J. 



B. Strong coupling regime 

In the strong coupling (SC) regime 3 ^> 1, equation (jSJ) shows that the lateral separation of counterions becomes 
larger than the GC length, thus indicating that the counterions tend to form a quasi-2D layer at the surface (Fig. Wp)- 
Such a layer is dominated by strong mutual repulsions between counterions as can be seen by considering the effective 
2D plasma parameter [56^ T = q 2 £B/a± ~ S 1 / 2 , which gives the ratio between Coulombic inter- particle repulsions 
and the thermal energy. For elevated S, Coulombic repulsions tend to freeze out lateral fluctuations of counterions on 
the surface, leading to strong correlations and a trend toward crystallization in the ionic structure [UEl- Individual 
counterions thus become isolated in relatively large correlation holes of size a± from which neighboring counterions are 
statistically depleted. The Wigner crystallization of the 2D one-component plasma is known to occur for T > T c ~ 125 



[561 ] . which corresponds to the range of coupling parameters S > S c ~ 3.1 x 10 4 [14 1. 

For S ^> 1, the PB description completely breaks down, nonetheless, one can obtain a simple analytical theory by 
employing a virial and 1/S expansion to the leading order, which is known as the strong coupling theory [l3|. The SC 
theory turns out to contain contributions that involve only single-particle interaction energies between counterions 
and the fixed macroion surface charges. For instance, the SC density profile of counterions at a single charged wall 
comes exclusively from the vertical degree of freedom, z, through which single isolated counterions are coupled to 
the wall with the interaction potential u/{ksT) = z/fi. Hence using the Boltzmann weight, one finds the (laterally 
averaged) density profile 

n sc (z) = n cxp(~z/fi), (9) 

where the prefactor (contact density) is again found from the electroneutrality condition to be n = 2-nt^a 2 . Unlike in 
the WC case, the SC density profile decays exponentially away from the charged wall. Moreover, the average distance 
of countcrions is finite and equal to the GC length, (z)sc = A*, reflecting again the quasi-2D structure of the layer. 

Formally, the single-particle SC theory is exact in the asymptotic limit of an infinitely large correlation hole size, 
aj_/ u — > oo, or simply S — > oo. However, its validity holds in a wider range of system parameters as is evident from 
comparison with computer simulations (l3.[T3|. For instance, for a finite coupling parameter 2, the SC density profile 
([9|) still holds approximately at distances z < a± , which yields the criterion 

z/n < VZ, (10) 
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a ± <ju<D 

D ^ ^ ^ ® d 

? @ q » a <a ® <a = = — = = = = = x 



(a) (b) 
FIG. 3: Schematic representation of the WC (a) and SC (b) interaction regimes for two like-charged walls. 



identifying the limits of applicability of the SC theory. At larger distances z > a±, multi-particle interactions become 
increasingly more important and the mean- field- like features eventually dominate for z»ai fl3l (3) ■ 

In brief, thus, one can identify two asymptotic regimes of weak and strong coupling based on the value of the 
electrostatic coupling parameter, where a charged system may be studied by means of two limiting theories, namely, 
the WC (mean-field) and the SC theory. 

In Table HI we present illustrative examples of the parameter values (surface charge density er s , counterion valency 
(?, GC length /x, and the coupling parameter S in water and at room temperature) from a few realistic weakly and 
strongly coupled systems. Note that a typical coupling strength of H ~ 10 2 (or larger) already falls within the SC 
regime and a value of S ~ 1 (or smaller) typically well inside the WC regime [13l - 4l6l | . 



V. INTERACTIONS BETWEEN LIKE-CHARGED SURFACES 



Macroions in solution are often like-charged and thus repel each other by their bare Coulomb interaction. The overall 
interaction is however different from this bare interaction due to the presence of countcrions, which can mediate both 
repulsive and attractive effective forces. Obviously, the counterion-mediated interactions depend strongly on the 
distribution of counterions around macroions. 

In order to demonstrate the underlying physical picture, we shall focus on the interaction between two identical 
planar like-charged walls of uniform surface charge density — cr s eo at separation D, where g-valent counterions fill 
only the space between the walls, Fig. [3] (the dielectric constant is assumed to be uniform in space). In this system, 
an extra length scale is set by the wall separation, D. Two limiting regimes of repulsion and attraction may be 
distinguished qualitatively by comparing D with other length scales of the system as follows. 



A. WC regime: Repulsion 



Let us first consider the WC limit, S -C 1, in the regime where the wall separation, D, is large compared with 
all other length scales (Fig. [3^,). In this case, counterions form a diffuse layer at each wall, but due to the large 
separation, the system is approximately decoupled into two nearly neutral sub-systems, each consisting of a charged 
wall and its counterion cloud. The effective pressure acting between the walls is dominated by the osmotic pressure of 
counterions across the mid-plane, which is positive and thus corresponds to an effective repulsion between the walls. 
The scaling behavior of this effective repulsion follows by noting that the mid-plane osmotic pressure can be estimated 
from the local counterion density, n m id, and by using the ideal-gas equation of state as f3P ~ n^^. Thus, according 
to Eq. (j8]), the interaction pressure is expected to decay as ~ D~ 2 . 

It turns out that the PB pressure obtained in the limit 3 — > coincides exactly with the mid-plane osmotic pressure 
of counterions as discussed above and may be expressed as @, H, HH 

^^) = A (11) 
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where A is obtained from A 1 / 2 tan[A 1 / 2 (D/2(ij\ = 1. From here one can obtain the large-separation D/fi 1 behavior 

9Ppb(D) „ 
2ir£ B a? ~ V D ) 



/3P PB P)^ m 2 (12) 



B. SC regime: Attraction 



Now let us consider the SC limit, S 3> 1, in the regime where D is smaller than the lateral spacing between 
counterion, di > D (Fig. Et>). Since counterions are highly separated from each other, the two opposite layers of 
counterions tend to form an inter-locking pattern at small separations. 

The system may be thought of as a collection of laterally frozen "correlation cells", each consisting of a single 
counterion sandwiched between two opposing sections of the walls with lateral size of about a±. Since aj_ 3> D, the 
effective pressure between the walls is dominated by the contribution from single counterions fluctuating in single 
correlation cells. The electrostatic energy of the system per cell is the sum of the bare interactions between the two 
surfaces with each other and with the single counterion, which-using the electroneutrality condition per cell and the 
fact that the wall separation is small-follows as u c \ cc / (ksT) ~ 2tt£b<J 2 D per unit area. This energetic contribution 
gives an attractive pressure as /3P c i cc ~ — 27t£bCs between the walls. On the other hand, the entropic contribution 
due to counterion confinement is of the order S C i ~ fee In D (per cell), which generates a repulsive component. The 
total pressure between strongly coupled walls is then obtained by combining these two effects as 

2n£ B a 2 ~ 1 + D ■ [U) 

This expression clearly predicts a closely packed bound state for the like-charged walls with an equilibrium surface- 
surface separation, D* , equal to twice the GC length, i.e., D* = 2fx. The like-charged walls therefore attract each 
other for D > D* and repel at smaller separations. 

The analytical expression (fT5]l is indeed an exact result for planar walls in the limit H — > oo [lU, [l4|. It turns 
out that in a system with finite coupling parameter, S, the asymptotic (S — > oo) SC results still hold approximately 
as long as the surface separation, £>, is smaller than the typical lateral distance between counterions, a±_, i.e., for 
D < a±. This condition in fact yields a simple and generic criterion identifying the regime where the SC attraction is 
expected to emerge between two like-charged macroions. It was originally suggested by Rouzina and Bloomfield (l7j 
and verified and generalized later using extensive analytical and numerical methods fl3l - fl6l |4ll I46l - l49l l5l| . For larger 
inter-surface separations, D > a±, the mean-field features become increasingly more imp ortant and the strength of 
attraction reduces. Eventually at very large D, the interaction becomes repulsive [l3|, Il4| . 



VI. COUNTERIONS WITH SALT 



The SC theory was so far designed exclusively for counterions-only systems, i.e., Coulomb fluids composed of only 
counterions in the absence of any salt ions (l3| . Though an approximation of this type can be used to describe situations 
where a large amount of polyvalent counterions dominate the system, it has to be amended in the general case in order 
to deal with the complexity of real systems that always contain some amounts of simple salt @ . An experimentally 
oft-encountered situation would be a system composed of fixed surface charges with polyvalent counterions bathed in 
a solution of univalent salt [57], [H[ . 

This situation leads to a difficult problem of asymmetric aqueous electrolytes where different components of the 
Coulomb fluid are differently coupled to local electrostatic fields [Hjj. Polyvalent counterions are coupled strongly, 
whereas univalent salt ions are coupled weakly. In this case no single approximation scheme that would treat all 
the charged components on the same level would be expected to work. Whereas the SC framework would certainly 
work for the polyvalent counterions, it would fail for the univalent salt. The converse is true for the WC framework. 
One is thus faced with a problem since no single approximation scheme appears to be valid in any range of coupling 
parameters. One can nevertheless build a theoretical framework that allows to selectively use different approximation 
schemes for different components of the asymmetric Coulomb fluid. This combined WC-SC approach appears to bring 
forth all the salient features of these asymmetric systems at high electrostatic couplings [4l[ . 



A. Functional integral formalism 



Our arguments until now were strictly intuitive. A formal theory can be developed exactly in terms of the functional 
integral representation of the classical partition function of the Coulomb fluid along the lines first introduced by 
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Edwards and Lcnard 

Assume first that the system is composed of charged macromolecules with fixed charge density po(r), mobile poly- 
valent counterions and an additional univalent salt. The total electrostatic interaction energy of a given configuration 
of the system can be written as 

W = i/" p(r)u(r,r')/5(r')drdr', (14) 

where v(r, r') is the Coulomb kernel given by v(r, r') = l/(47reeo|r — r'|), and p(r) is the total charge density 

P (r) = Po (r) + £ qe 5(R%- r) + ]T e S(R+- r) - £ e 5(Rr- r), (15) 

i i i 

where R^, and Rr are the positions of the polyvalent counterions, univalent cations (salt counterions) and 
univalent anions (salt coions), respectively. 3 Assuming again that the system is composed of two apposed planar 
surfaces at z = ±D/2 with the surface charge density — ct s , we have 

p (r) - -a s e [5{z - D/2) + 5{z + D/2)] . (16) 

The salt ions are assumed to be present in all regions in space, whereas the counterions are assumed to be present 
only in the inter-surface region \z\ < D/2 as specified by the geometry "blip" function Sl(r) (Section pv Al) ). 

We then follow the standard procedure by introducing a fluctuating local potential, 0, via the Hubbard-Stratonovich 
transformation, which leads to the following exact functional integral representation for the grand-canonical partition 
function [H, [l]| 



V</>e- pH[ *\ (17) 



where the field-functional Hamiltonian reads 



H \<t>\ = \J 0(r>- 1 (r,r / )0(r / )drdr / + i J Po(r)<^(r)dr (18) 



P J P J P 

and A c and A± represent the fugacities of polyvalent counterions and salt ions and w _1 (r, r') = — e£qV 2 5(y — r') is the 
inverse Coulomb kernel. The special case of counterions-only system, as analyzed in the previous section, is recovered 
by setting A± = 0. 

We shall assume that salt ions are in equilibrium with a bulk reservoir containing equal concentration rib of both 
positive and negative ions, which implies A + = A_ = rib. One can thus introduce the Debye-Hiickel (DH) screening 
parameter k (inverse "screening length") as k 2 = 87t£b"-6- The polyvalent counterions shall be treated here within the 
canonical ensemble 4 , i.e., their number in the slit is assumed to be fixed and equal to N. The number of counterions 
may be expressed via the dimensionlcss parameter 

r, = Nq/(2a s S), (19) 

where S is the area of the interacting surfaces. The case rj = represents a system with salt only, and rj = 1 is the 
case where the total charge due to counterions exactly compensates the surface charge. Note that r\ can take any 
non-negative value when salt ions are present. This is because salt ions turn the long-range Coulomb potential into 
a short-range DH potential (see below) and can thus ensure the electroneutrality condition themselves. 

B. Dressed counterions 

Assuming that the system is highly asymmetric g> 1, one can formulate an approximate theory in order to evaluate 
the partition function (| 1 T[) analytically by acknowledging the fact that the polyvalent counterions are strongly coupled 



3 We may refer to the q-valency (polyvalent) counterions simply as "counterions" 

4 see Ref. [4lll for a grand-canonical description of polyvalent counterions. 
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while the simple salt ions are weakly coupled to the fluctuating electrostatic fields. This leads to a mixed WC-SC 
evaluation of the partition function [4l| . 

The salt terms (the last two terms) in Eq. (|18j) can be combined into cos /3eo0(r) and in a highly asymmetric 
system this can be expanded up to the quadratic order in the fluctuating potential. Thus up to an irrelevant constant 
we remain with an effective field Hamiltonian of the form 

H eS [^} = l J 0(rH^(r,r')0(r')drdr'+i J p (v)c/>(r) dr 

_^ [ c - i/3 « en< ^(r)dr. (20) 
P J 

This procedure therefore yields an effective Hamiltonian for a "counterions-only" system but with the proviso that 
the inverse Coulomb kernel is replaced by the standard inverse DH kernel 

^(r,r')--££ (V 2 - K 2 )<5(r-r'), with « DH (r,r') = (21) 

We have thus effectively integrated out the salt degrees of freedom leading to a renormalized interaction potential 
between all the remaining charge species of the screened DH form. One can thus drop any reference to explicit salt 
ions and infer the thermodynamic properties of the original system by analyzing it as a system composed of dressed 
counterions and fixed external charges interacting via a screened DH pair potential. In the SC limit for the polyvalent 
counterions we term this approximation scheme as the SC dressed counterion theory. Our SC analysis thus proceeds 
in the same way as for the counterions-only systems [l3| except that the interactions between the charges arc now of 
a dressed form. 

We note that any Bjerrum pairing l60l | or even electrostatic collapse of the salt or formation of salt-counterion 
complexes [6l| is beyond the framework developed here. 

C. WC dressed counterion theory 

We again focus on a system composed of two plane-parallel surfaces defined via Eq. ([T6| . In the WC limit (for 
both the counterions as well as the salt ions), the functional integral derived in the previous section is dominated 
by the contribution from the saddle-point solution 0sp- This subsequently leads to the mean- field equation for the 
real- valued mean- field potential ip = i'/'SP, i-e., 

- ee (V 2 ^ - k 2 ^) = Po (r) + qe A c Q(r) e ~^, (22) 

which, in rescaled units w — f3qeoip and by virtue of the lateral symmetry for planar surfaces w ~ w(z), may be 
written as 

w" = k 2 w - Ce~ w \z\ < D/2. (23) 

The constant C can be evaluated when one stipulates the fixed amount of counterions. Outside the slit \z\ > -D/2, 
the mean-field equation has the standard DH form w" = k 2 w, which yields w(z) = wq exp(±«z). 

The interaction pressure, P, between the bounding surfaces is given by the difference of the ion concentrations at 
the mid-plane (z = 0), where the mean electric field vanishes, and the bulk concentration, i.e., (3P — n + (0) +?t._(0) + 
n c (0) — 2rih, which leads to the dimensionlcss expression 

|^4 = i(^) 2 w 2 (0) + ice-^). (24) 

As evident from the above equation, the pressure can never be negative and the effective interaction is thus always 
repulsive within this type of mean-field approach (62l | . The canonical mean-field equation (|23p can be solved numerically 
and the results can be used to evaluate the pressure (|24[) as a function of the inter-surface separation (Fig. [4^,, dashed 
line). 

D. SC dressed counterion theory 

The analysis of the dressed counterions system in the SC limit is very similar to the traditional SC approach in the 
case of counterions only (l3j . We proceed by expanding the grand-canonical partition function associated with the 
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dressed counterion approximation, Eq. (|20[) , to the first order in counterion fhgacity, A c . We then perform an inverse 
Legendre transformation 5 in order to obtain the canonical SC free energy |l3l. l46j. We thus find 

F N = U - Nk B T In J c~ Pu( - r) dr, (25) 
where the first term is the screened interaction energy of fixed charges 

Uo = \\ p (r)«DH(r > r')po(r')drdr' I (26) 

and the term in the exponent is the single-particle interaction energy of the dressed counterions with fixed macroion 
charges 

u(r) = qe J p {r')v DH {r, r') dr'. (27) 
The SC attraction between like- char ged macroions stems from the second term in Eq. (|25|) . which contains the 



the bU attraction between like-char ged 
counterion-induced effects [HI, HH, El. 

For the planar system under consideration (Eq. (|16[) ). the above quantities may be evaluated explicitly (4l| and we 
find the SC free energy as 

PFn/N = c- KD - In 1(D), (28) 

where we have introduced 

fD/2 2 

1(D) = / cxpf — c- KD/2 cosh Kz)dz. (29) 

J \Kfl J 

The first term in Eq. (|28[) corresponds to the usual salt-mediated repulsive DH interaction between the two surfaces, 
and the second one is the contribution of counterions, which is proportional to 77, Eq. (|19|) . on the single-particle SC 
level. 

The dimensionless pressure acting between the surfaces can be obtained from the free energy via the standard 
thermodynamic relation P = — 8(Tm /S)/dD, thus yielding 

where the prime denotes the derivative with respect to the argument. 

The analytical SC pressure ([50)) is shown in Fig . @Jt as a function of the inter-surface distance (solid curve) along 
with the results from Monte-Carlo simulations [41[ (symbols) for a few different screening parameters. We also show 
in the figure the WC pressure from Eq. ([24)) (dashed curve). First note that the simulation results are bracketed 
by the two limiting analytical theories of WC and SC within the dressed counterion scheme and thus agree with the 
general feature obtained before EMU that the WC and SC limits in fact establish the upper and lower bounds for 
the interaction pressure between charged surfaces. 

For both small and larg e sep arations the interaction pressure becomes positive (repulsive), whereas for a sufficiently 
large coupling parameter [4l| . an effective inter-surface attraction can emerge at intermediate separations between 
the walls. For kD <C 1, the SC pressure reduces to 

__ * — + (1 - 2 V ), (31) 

which to the leading order corresponds to the ideal-gas osmotic pressure of counterions squeezed between the two 
surfaces. For kD >• 1, the pressure behaves as 

2n£ B ^ ~ D (kD) 2 ' 1 ' 



5 This is achieved by mapping the fugacity to the number of counterions, iV, via the relation A c 91n Z/dA c = N. 
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D/jU D/\i 

(a) (b) 

FIG. 4: a) Rescaled interaction pressure P = /]P(D)/(2iy£-b(Ts) between two like-charged surfaces as a function of the rescaled 
inter-surface distance, b) Van der Waals type iso-inverse screening length curves are shown along with the corresponding 
Maxwell construction. The inverse screening length is varied in the range from k/i = 0.4 to k/i = 0.56 in intervals of 0.027 
(from bottom to top) . 



which indicates that at large separations the counterions again behave as an ideal gas on the leading order as all 
electrostatic interactions arc effectively screened out and hence only the repulsive osmotic contribution remains. 

The interaction pressure in the canonical ensemble thus always possesses repulsive branches at small and large 
separations and can show non-monotonic behavior in between. In fact, for certain values of the parameters the 
interaction pressure shows a van der Waals-like loop which could suggest a coexistence regime between two different 
"phases" . This loop is obtained for certain iso-ionic strength curves. From a thermodynamic perspective one thus has 
a coexistence between a dense phase, identified with a small inter-surface separation, at equilibrium with an expanded 
phase with a larger inter-surface separation. Such van der Waals-like coexistence between interacting charged surfaces 
has been seen in other contexts before (63j and can be demonstrated by means of a Maxwell construction analysis as 
shown in Fig. @J). 

In an experiment such as osmotic stress experiments and surface force experiments @ one can only probe stable 
equilibrium states of the system implying interaction pressure vs. separation curves that are in agreement with the 
appropriate Maxwell construction. The binodal or the coexistence curve, which delimits the region in the pressure- 
separation plots where a Maxwell construction is feasible (red curve in Fig. Hp), ends at a critical point corresponding 
to a critical amount of salt above which the interaction pressure remains purely repulsive. For the case with r\ = 1, 
we find the critical point as (P c = 0.092, D c /^l = 6.14, /t c /x = 0.546). 

It is interesting to note that this type of interaction pressure equilibria corresponding to abrupt transitions from one 
equilibrium separation to a different one have been observed in experiments with strongly charged macromolecules in 
the presence of polyvalent counterions and univalent salt. A typical example would be the osmotic stress experiments 
on DNA in the presence of trivalent CoHex counterions and 0.25M NaCl salt, that show abrupt transitions in osmotic 
pressure for intermediate CoHex 3+ concentrations from one repulsive osmotic pressure branch to another one 57j |. 
Similar features are discerned even for a divalent counterion Mn 2+ at various concentrations or temperatures [58 1. 

The agreement between the SC dressed counterion theory and simulations becomes better as the coupling parameter 
S becomes larger. The agreement is also better for a smaller fraction of counterions, ?y, in the slit. Using a similar 
argument as in the counterions-only case in Section IV Bl [13l-[l6j , we find that the theory is expected to hold at small 
separations given by D//i <C y/'B./rj. Thus, for r\ < 1 (i.e., when the amount of the bare charge due counterions is 
less than the bare fixed charge on the macroions), the SC dressed counterion theory is expected to hold in a wider 
range of separations as compared with the original counterions-only SC theory At very large separations, where 
most of the electrostatics is screened out, the interaction between counterions becomes negligible and the SC theory 
of dressed counterions retains its validity again. This result is a consequence of the dressed counterion theory and is 
not obtained in the standard SC theory with counterions only. Thus, the SC dressed counterion theory captures the 
physics both at large and small separations but would require improvements at intermediate separations. 

It should be noted that the validity of the DH-type linearization that we have used to derive the dressed interaction 
potentials is also limited by stipulating that the dimcnsionless DH potential itself is always small enough. This leads 
to the condition that n 3> 27r€BC s , or k/i > 1/g, which turns out to cover a whole range of realistic parameter values 
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VII. COUNTERIONS BETWEEN RANDOMLY CHARGED SURFACES 



The assumption of homogeneity of surface charges is in general quite severe and there are well known cases where 
this assumption is not realistic at all. Random polyelectrolytes and polyampholytes present one such case [64T[66| . 
There the sequence of charges can be distributed along the polymer backbone in a disordered manner where the 
disorder distribution may be of a quenched type. The Coulomb (self-)interactions of such polyelectrolytes are distinct 
and different from homogeneously charged polymers. 

A case even closer to the present line of reasoning are investigations of interactions between solid surfaces in the 
presence of charged surfactants. The aggregation of surfactants on crystalline hydrophobic substrates in aqueous solu- 
tions can sometimes show structures consistent with highly inhomogeneous and disordered surface charge distributions 
[67j . Similar interfacial structures are seen for interacting hydrophilic mica surfaces in the presence of cetyl-trimethyl- 
amonium bromide (CTAB) or other surfactant-coated surfaces [68| . The surfaces appear to be covered by a mosaic 
of positively and negatively charged regions and experience a strong, long-ranged attraction, which is comparable in 
magnitude to that between hydrophobic surfaces, and is orders of magnitude larger than the expected Lifshitz-van der 
Waals forces [(38| . The patterning of interacting surfaces described above is highly disordered, depends on the method 
of preparation and has basic implications also for the forces that act between other types of hydrophilic surfaces with 
mixed charges. 

It thus seem appropriate to investigate the effect of quenched disordered charge distribution on the interactions 
between macroions in ionic solutions. 6 As a particular case, we shall again focus on the effective interaction between 
two randomly charged planar surfaces across a one-component Coulomb fluid [42j . 



A. General formalism: The replica method 

The partition function of a Coulomb fluid in the presence of an external fixed charge distribution pa(v), can be 
again written in the form of a functional integral over the fluctuating electrostatic field <f>(r) as given in Eq. (|17[) . 
However, the fixed charges are now assumed to be randomly distributed on macromolecular surfaces. Thus po(r) is 
represented by a probability distribution, which is assumed to be Gaussian with no spatial correlations, i.e. 

V[ Po (r)} = const, x A I dr (Po(r) - PoMf (33) 

where po( r ) is the mean value and <?(r) the width or variance of the charge disorder distribution. For clarity, we shall 
also focus on the counterions-only case by setting A± = in Eq. (fT8|) . 

The average over quenched charge disorder is now obtained by applying the standard Edwards- Anderson replica 
ansatz (gl, in the form 

~2™ _ i 

T = -k B T\nZ = -k B T lim , (34) 

n->0 n 



where the disorder average is defined as (...) = J Vp (• • • )V\po{r)]. 

The Gaussian integrals involved in Eq. p4j) can be evaluated straightforwardly and the final form of the replicated 
partition function follows as [13] 



e -/J5[*.(r)] j (35) 



with 



%(r)] = ^ / 0a(r)P a6 (r,r')^(r)drdr'+i / p (r) £ o (r) dr 

a.b ^ J a 

~ y*r!(r)^e- i/W »Mdr, (36) 



See Refs. [44l, 145H for an analysis of the effects due to annealed and partially annealed charge disorder. 
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where a, b = 1, . . . , n are the replica indices and 

V ab (r,r') =u- 1 (r 1 r')6 ab + f3g(r)5{r-r'). (37) 

The expression ([35]) together with Eq. (|34|) represents the starting formulation for the free energy in the presence 
of quenched charge disorder. This free energy can only be evaluated approximately. Thus, in order to proceed we 
shall combine the methods developed for the one-component (counterions-only) Coulomb fluid without disorder (l3j 
and modify them in order to incorporate appropriately the disorder effects. We shall start with the WC limit giving 
rise to the corresponding mean-field theory and then proceed to the SC limit. 



B. Disorder effects in the WC regime 

In the WC limit 5 — > 0, one may proceed by employing a saddle-point analysis of the functional integral Q35p . just 
as in the case with no disorder in Section IVI CI It is easy to show that the real- valued mean-field replica potential ip a 
is governed by the following equation, 

n 

-ee V 2 V; a (r)+/? 5 (r)5>b(r) = Po(r) + ge A c fi(r) e"^-. (38) 

6=1 

In the replica formulation we have to take the limit n — > 0, which furthermore implies that linin^o ^Pbi?) — ^ 0. It 
is thus evident that in the limit n — > 0, the contributions from the disorder vanish and, because the index a becomes 
irrelevant, one recovers the standard PB equation 0. Therefore, the quenched charge disorder effects completely 
vanish in the WC limit 0, |Hj]. 

The above result is a consequence of the mean-field approximation and holds in the limit 3 — > even if the system 
is generalized to contain additional ionic species or dielectric discontinuities at the bounding surfaces. The quenched 
charge disorder however turns out to play a significant role in dielectrically inhomogeneous systems when electrostatic 
field fluctuations are taken into account. It can be shown to lead to an additional attractive or repulsive contribution 
to the total free energy (depending on the dielectric mismatch and the salt screening in the system) even when the 
surfaces are assumed to be net-neutral [H, [44| . 



C. Disorder effects in the SC regime 

The partition function (|35p can be calculated in the SC limit via a virial expansion up to the first nontrivial leading 
order in powers of the fugacity as noted before. The canonical SC free energy of the system then follows from Eq. 
(|34|) by using a standard Legendre transform [42j as 

~ FN= lf ^o(i-)w(r,r')p (r')drdr' + iTr.g(r)t;(r,r') 

-Nk B T In / fi(R) c-P u{¥t) dR. (39) 



The first term in Eq. (|39[) is nothing but the direct Coulomb interaction between the mean charge densities po(r) of 
the fixed charged surfaces (macroions). The second term is an additive contribution from the charge disorder, which 
becomes important only in dielectrically inhomogeneous systems [43l |44j and will be irrelevant in the present study. 
The third term, however, embodies the disorder effects on the SC level in the presence of a Coulomb fluid. It depends 
on the single-particle interaction potential 

u(R) = qe Q J w(r', R)po(r') dr' - ^(qe Q ) 2 J g(v')v 2 (r', R) dr', (40) 

where the second term comes from the disorder variance g(r) and exhibits a non-trivial dependence on the Coulomb 
kernel v(r,r') = l/(47ree |r — r'|). 

Assuming again that our system is composed of two planar surfaces located at z = ±D/2 with statistically identical 
random charge distributions, we write the mean density and variance of the disordered surface charge as 

p (v) = -a s e [S(z-D/2) + S(z + D/2)}, g(r) = ge\ [5(z - D/2) + 5{z + D/2)} . (41) 
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(a) (b) 

FIG. 5: Left: Rescaled SC free energy, Eq. (|42[) . of two charged walls bearing quenched charge disorder as a function of the 
rescaled inter-surface distance D/p, for \ — 0, 1 and 2 (bottom line, middle line, top line). Right: Rescaled equilibrium distance 
D* I \i as a function of the disorder coupling parameter \. 



The electroneutrality again stipulates that 2a s S = Nq. The geometry function f2(R) is the same as before. All the 
terms in the expression for the free energy, Eq. (|39p , can be computed explicitly. At the end we obtain a surprisingly 
simple expression 

^ = - + ( X -l)lnD. (42) 

Here we have introduced the dimcnsionlcss disorder coupling parameter 

X = 2irq 2 £ 2 B g, (43) 

which is very similar to the electrostatic coupling parameter, S, in Eq. (j4]), except that it is defined based on the 
disorder variance g and scales with the countcrion valency as q 2 instead of q 3 . The free energy (|42[) is plotted in Fig. 
[5ji for different values of the disorder coupling parameter. Note that the disorder leads to a long-range attractive 
contribution, which is additive in the SC free energy and has a logarithmic dependence on the separation, i.e., ;\Tn ZJ. 
It thus appears that the quenched charge disorder and the counterions confinement entropy, i.e., the — In D term in 
Eq. (|4"2"j) . in some sense counteract one another. 

Evaluating the interaction pressure from the free energy, Eq. (|42|) . we find P(D) = Psc{D) + -Pdisorder(-D), where 
the first term is the standard SC pressure [Hj], Eq. ([T5)|. and the second term is the additive contribution from the 
disorder 

^disorder (J) _ /2/A 

2ne B o-z X \Dj- 1 > 

We can then derive the equilibrium distance D* between the two surfaces, corresponding to zero interaction pressure, 
as 

D* = 2(l- X )/i. (45) 

In the undisordered case, x = 0, this reduces to the known result D* = 2/i [T3|, which corresponds to a stable bound 
state for the two surfaces at a separation equal to twice the GC length. However, as x i s increased, the equilibrium 
separation decreases and vanishes at the critical value Xc = 1 and remains at zero thereafter. This behavior has 
all the features of a second-order, qucnchcd-disordcr-induccd collapse transition with an unusual value of the critical 
exponent (see Fig. [5}j). Note also that for X = 1 the interaction pressure between the surfaces is obviously constant 
in the whole range of separations D right down to zero as the counterions confinement entropy is completely wiped 
out by the charge disorder contribution. 
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VIII. LESSONS 

The two limiting laws, i.e., the WC and the SC limits for Coulomb fluids that we explored above, are valid in disjoint 
regions of the parameter space. While the WC limit is valid for sufficiently small macroion surface charge densities, 
low countcrion valencies, high medium dielectric constant and/or high temperatures, the SC limit becomes valid for 
respectively opposite parameter values. The two together bracket the region of all possible behaviors of Coulomb 
fluids confined between charged boundaries, a view that was completely corroborated by extensive simulation studies. 

The parameter space in between these limitin g va lues can be analyzed by appr oximate methods [724474^ but is most 
often accessible solely via computer simulations jl3l4lrl [ill I46H491 [Bll. l72l |744476| . Exact solutions for the whole range 
of coupling par ameters are unfortunately available only in one dimension [77] . The WC-SC paradigm has been tested 
extensively jili - [ltl l4ll . [46M49L IFll [721474 [l3] and fits computer simulations quantitatively correctly in the respective 
regimes of validity, thus providing a unifying conceptual framework of the behavior of Coulomb fluids. 

Though we have shown that in the important limit of SC the much cherished and widespread PB approach docs 
not work, formally its applicability can be systematically extended by perturbative corrections in the local potential 
fluctuations and correlations [1343 [EH [EH l78l48l| along the lines of the standard approach used in the mean-field 
context [82j ■ This kind of fix is nevertheless severely limited since the perturbative expansion is only weakly convergent 
[llHlH and higher-order corrections beyond the first-loop Gaussian term are very complicated and difficult to carry 
through [8(| l8]J|. S uch perturbative corrections offer in effect only a relatively insignificant improvement over the PB 
approximation [I3l4l6[ l5l| and can not predict phenomena such as like-charge attraction [75| . 

The emerging world of Coulomb interactions reviewed above is indeed fascinating. While the counterion-mediated 
electrostatic interactions between equally charged surfaces are always repulsive on the WC level, the SC regime offers 
a much richer framework with plethora of new phenomena. Interactions between equally charged macromolccular 
surfaces can be attractive for strongly coupled counterions, or non-monotonic-showing repulsion at small and large as 
well as attraction at intermediate separations-for strongly coupled counterions in the presence of weakly coupled simple 
salt. While charge disorder on macromolecular surfaces has no effect on the WC level, it can quite unexpectedly lead 
to strong electrostatic attractions between randomly charged surfaces on the SC level. In view of these developments, 
the commonly held pop culture wisdom that likes repel and oppositcs attract should thus be substantially amended! 
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